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TECHNICAL MEMORANDUM X-64535 


GENERAL RELATIVISTIC PRECESSION OF A GYROSCOPE 

IN AN INCLINED ORBIT 

SUMMARY 


The general relativistic precession of a gyroscope in an inclined orbit 
is calculated. The magnitude of the precession is found to be proportional to 
the cosine of the inclination angle, apart from periodic terms. It is shown 
that the geodetic precession cannot be separated from the motional precession 
no matter what initial orientation is chosen for the gyroscope. 


INTRODUCTION 


A perfectly spherical gyroscope in orbit around the earth will precess 
as a result of two relativistic effects. The first and larger effect results from 
the motion of the gyroscope along a geodesic in a 4-space which is not flat; 
this is called the geodetic precession. The second effect is associated with 
the rotation of the earth and is called the motional or Lense-Thirring preces- 
sion. A polar orbit was originally chosen for the Stanford Gyroscope Relativity 
Experiment because this orbit will allow these two effects to be separated. 
Recently, it has been suggested that a preliminary test flight of the Gyro 
Experiment be performed on Skylab II, for which an inclined orbit is planned. 

In this paper we calculate the expected gyroscopic precession for such an 
inclined orbit and show that one cannot separate the geodetic and motional 
precession. 


THE GEODETIC PRECESSION 


First, we will consider the geodetic precession, which is generally 
about two orders of magnitude larger than the motional precession. To obtain 
this precession we must solve Schiff's equation [1] for the orbit under con- 
sideration. This equation will be written as 


S x L 




dS 

dt 


(1) 





1 GM. 

9 ^ 2^3 

v X 


(r x v) 


( 2 ) 


which is a first order equation in the gyroscope spin vector S. 2GM/e 2 is the 
Schwarzschild radius of the earth, and~r and v are the orbital position and 
velocity vectors. L has been defined to be antiparallel to the orbital angular 
momentum vector. 

In the simple case of a circular orbit for which L is constant, the 
vector S simply precesses about L with constant angular velocity c o - | L | 
and S • L remains constant. For the more general case in which L is not 
constant, one must obtain T and v, thus L, as functions of time from an 
integration of the orbit equations and then obtain S as a function of time by 
integrating equation (1). 

For a circular inclined orbit about the earth, the orbital plane will 
not remain fixed in space. Because the earth is not perfectly spherical, the 
orbital plane will precess about the earth's axis in such a way that the orbit 
inclination angle and radius remain nearly the same. So to good approxima- 
tion, one can consider that the orbital angular momentum, and thus L, will 
remain constant in magnitude and precess about the earth's axis with an angular 
velocity fi. Classical perturbation theory [2, equation 11. 15.6] gives fi as 


n 



(cos i) 


C 0 

o 


(3) 


where i is the orbital inclination, is the orbital angular velocity, and J 2 is 
the earth's mass quadrupole moment. L can then be written in the form 


L = Lm sin fit 
x 


L = Li cos fit 

y 



(4) 


where and L 2 are positive constants. The z-axis has been chosen to 
correspond to the earth's axis while the x-y plane coincides with the equatorial 
plane. L has been given a negative z-component so that the orbital angular 
momentum vector has a positive z-component. Initially L is in the y-z plane. 
It precesses in a direction consistent with the fact that the nodes regress for 
an oblate body. 
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When L is given by equation (4) , the following analytic solution to 
equation (1) can be obtained, 

S = -C 1 sin p sinfit - C 2 [sin (cot + i p) cos fit - cos (cot + tp) sinfitcos p] 

S y = -Ci sin p cosfit + C 2 [sin (cot + ip) sinfit + cos (cot + ip) cosfitcos p] 

S = C 2 sin p cos (cot + p) + Ci cos p , (5) 

z 

where 

co = ^L 2 + fi 2 + 2L 2 fiJ 1/2 

sin p = L^/ co , 

cos p = (L 2 + fi)/c o , 

and the constants Ci, C 2 , and p are determined by the initial condition on S. 

We will now demonstrate explicitly that this is a solution to equation (1) . By 
differentiating equations (5), we obtain 

dS 
x 

dt 
dS 

y 

dt 
dS 

z 
dt 

-£2Ci sin p cos fit + C 2 [(fi - co cos p) sin (cot + i p) sin fit 

- (co -'fi cos p) cos (cot + p) cos fit] 

fiCi sin p sin fit + C 2 [ (fi - co cos p) sin (cot + p) cos fit 

+ (co - fi cos p) cos (cot + p) sin fit] 

- coC 2 sin p sin (cot + p) 


dS 

dt 


-[* 


+ Lo 


( 6 ) 
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Using equations (4) , we have 


(S x L) = 


‘s L 

- S L 


S [ — L 2 ] - S [Li cos fit] 

y z 

y z 

z y 


S L 

- S L 


S [Li sin fit] - S [-La] 

z X 

X z 


z X 

S L 
X y 

- S L 

y x 


S [Li cos fit] - S [Li sin fit] 
x y 1 


and now substituting from equation (5), this becomes 


Cj [L 2 sin 0 - Li cos <p] cos fit - C 2 L 2 sin (cot + ip) sin fit 
-C 2 [L 2 cos <p + Li sin cpj cos (cot + ip) cos fit 

-Ci [L 2 sin 0 - Li cos <p] sin fit - C 2 L 2 sin (cot + ip) cos fit 

+C 2 [L 2 cos (p + Li sin <p] cos (cot + ip) sin fit 

-LiC 2 sin (cot + ip) [cos 2 fit + sin 2 fit] - LiC 2 cos <p cos (cot + ip) 

[-sin fit cos fit + cos fit sin fit] 


Now from equations (6), we have 

L 2 sin <p - Li cos ^ = L 2 - Li = = -fi sin <p 


L 2 cos <p + Li sin <p = L 2 ( ^‘ 8 + ^ 


CO 


+ Li(^-) 

\ CO / 


(Li 2 + L 2 ) 


CO 


^ +fi(^ 

CO 


(co 2 - fi 2 - 2L,fi 


CO 


^ + fif^ 

CO 


CO - fi 


L 2 + fi 

CO 


= CO - fi COS (p 


-fi + (L 2 + fi) = - (fi - CO cos <p ) 
co sin <p , 


so that a comparison of the expressions for dS/dt and S x L shows that they 
are equal. Hence, equations (5) have been shown to be a solution of 
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equation (1) when L is given by equation (4). It can also be shown, using 
equations (5) that 

s’ • "s - S 2 + S 2 + S 2 - Cl 2 + c, 2 
x y z 

so that S remains constant in magnitude, as it must. 

This solution to equation (1) was obtained by transforming to a 
rotating coordinate system in which L is at rest. In this system, ~S precesses 
about L and one can write a solution to equation (1) in the rotating system 
immediately. When this solution is transformed to the nonrotating system, 
equations (5) are obtained. The problem is mathematically identical to the 
motion of a magnetic dipole in a magnetic field with a constant z-component 
and a rotating component in the x-y plane. This configuration is often 
considered in NMR work. 


Using standard identities for the products of two trigonometric 
functions, equations (5) can be written in the form 

S = -Cj sin <p sin fit - C 2 + - ■ s i n [ ( W _ fi)t + ip] 

-C 2 — ^° S ^ sin [ (co + fi)t + ip] 

Li 

S = -Ci sin cp cos fit + C 2 ■ ■ 1 cos [ (co - fi)t + ip] 


-C 2 ii ?2- s -.&) .- cos [ (k, + n)t + ip] 

Z 1 


S = Cj cos <f> + C 2 sin <fi cos (cot + ip) 
z 

so that S consists of terms with frequency co - fi, fi, co, and co + fi. We can 
simplify this expression by using the fact that L « fi; that is 


T 3GM | r x v | 

— o ~ ™~*> 

2c z r r 


2c R V r 


3GM /R 

— I co 


= [1. 05 x 10“ 9 ] 


R 

' co 

r / o 


6.6 x 10 


, 9 /R\ rad 

{r J 


rev 
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where R is the earth's radius. Now since J 2 ~ iO -3 .U , we have from 
equation (3) that £2 ~ 10" 3 co ~ 5 degrees/day for a near-earth orbit so that 

O 


ft 


1CT 9 

10" 3 


CO 

0 


CO 

o 


NT 6 


Expanding the first of equations (6) to first order in L/S2, we have 

r2l% 


CO 




< 2L™ 1/ 

1 + -=* + — 


n 


£2 


n 


, L 2 
1 + -~ 
n 


so that 


co — £2 — L 2 ( V ) 

and 

sin cp c* Lj/O , (8) 

all valid to first order in L/£2. By expanding co to second order in L/£2 and 
using the expression for cos cp in equations (6) , we find 



so that (1 - cos cp) is second order in (L/£2). 

If we let the initial value of S be S°, then at time t = 0, equations (5) 

give 


S = -Co sin ip 
x d r 


S y = -C* sin cp + C 2 cos ip cos cp 


( 10 ) 


S = Ci cos <p + C 2 cos ip sin cp 
z 


These three equations can be solved for Cj , C 2 and ip to give 


ip = tan 1 I -S / 
I x 


r , o , _ o . 

S cos cp + S sm cp 

y z 


6 



o 


( 11 ) 


n. 


S cos (b - S sin d> 

z y 


S X°) + ( S y°) cos2 ^ + ( S Z°) sin2 ^ 


+ 2 ( s y j( s z ° ) cos sin 0 


V 2 


and a substitution of equations (10) into equations (11) verifies that this is a 

solution. So that for any initial S°, a solution can be constructed. At any 
time t, the angle 0 through which the spin vector S has precessed can be 
obtained from the equation 



where S is given by equations (5) and S° by equations (10). This is not a 
very convenient way to compute 0. Since 0 is very small, it is easier to use 
the following method. Define the vector B by the equation 


(12) 

V2 

IB x ~S ° | _ [(BxS°) • (Bx~S°)] 

S 2 S 2 

Using a standard vector identity this becomes 
sin© = [S 2 B 2 - (S° • B) 2 ] 1/2 /s 2 

since 

(S°) 2 = (S) 2 . 

Now 

S 2 = (S° + B) ' (S° + B) = S 2 + 2 B • "s ° + B 2 


S = S + B 


Then, we have 


sin 0 


|S x S 
S 2 
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so that 


(13) 


(S • B) - -B 2 /2 


sin Q 


S 2 


» 2 B 4 

B 

i_ / jb\ 2 

S 2 B 2 + — 

__ 


4 

S 

+ 4 \S/ _ 


V 2 


(14) 


A general formula for sin 0 can be calculated using equations (5) , (10) , (12) , 
and (14) but the result is quite complicated. So we will restrict ourselves 

to three specific choices of S° and work out explicit expressions for 9 in 
these cases only. 


1. S y ° = S, S x ° = 0, S° = 0; i. e. , S is initially perpendicular to 

the earth's axis and in the plane formed initially by L and the earth's axis. 
From equations (10) we have 

tp = 0 , Ci = -S sin 0 , C 2 = S cos (p 


so that equations (5) become 


S 


x 


S 

y 


s 

z 


S 


s 


s 


sin 2 0 sin S2t - cos 0 ^ + ~ 3 sin (co - fi)t 

, (1 - cos 0) . . _ . , ~ 

- cos 0 sin (co + fi) t 

■ 2 . , (1 + cos <b) , 

sim 0 cos fit + cos 0 ^ cos (co - fi)t 

, (1 - COS 0) . 

- COS 0 T — cos (cu + fi)t 

z 

[cos 0 sin 0 (cos cot - 1)] 


Now from equations (8) and (9) we note that sin 2 0 and (1 - cos 0) are 
second order in L/fi. Also, we are interested in times of the order of 1 year 
so that from equation (7) 


(co - fi)t =* L 2 t « 1 


(15) 
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since L ~ 3. 5 x 10 -5 rad/year. We can then expand the sin function to obtain 


S, <* - S [ L 2 t] 

X 


(16) 


so that S is linear in t to first order in L/12. Also 

x 


S =* S 
z 


ill 

n 


[cos cot - 1] 


(17) 


from equations (8) and (9), so S oscillates with frequency co ~ 12 to first 

z 

order. The precession angle 0 can be calculated as follows: Using equations 
(12) and (13) we have 


B 

S 


„ -?o, 


-BS 

-2 (B • S t ) 

' l 

./o 

y y 

S 2 

— \J Z 

s 2 


\[2 


. 2 , , (1 + COS 0) , 

sin 0 cos 12t - cos 0 — — ^ cos (co - 12) t 


, (1 - COS 0) . ^ , , 

+ COS 0 C0S + 


If we now use equations (8) , (9) , and (15) to expand this to second order in 
L/12 and (co - 12 )t, we find that B/S « 1 and 


0 = I = [Ljtl! + (2 


3 1 

— - 2 cos fit + — cos (co + 12) t 
z z 


V 2 


(18) 


so 0 is not strictly linear in t but only becomes so for t >> 7^ ( 7- . . 

B 2 \12 / B 2 \27T 

=* 10 days, where t ~ 60 days is the precession period of 12. 

2. sj 3 = S, S ° = 0, sj 3 = 0; i. e. , S is initially perpendicular to 
the earth's axis and to the initial value of L. From equations (10) we have 

0 = -71-/2, C t = 0, C 2 = S 


(1 + COS 0) , (1 - COS 0) . . , 

^ cos (co - 12)t + — — cos (co + 12)t 
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(1 + cos <p) . . . (1 - cos d>) . , 

sm (co - O) t - — y sin. (co + fi)t 


S = S sin d> sin cot 
z 


Again we see that, since (1 - cos <p) is second order in L/fi, we have after 
expanding sin (to - fi)t to first order that 


S y - S [L 2 t] 


(19) 


and S has an oscillatory behavior with frequency co ~ fi. The precession 
z 

angle is calculated in the same way as is case 1: 

JVi 

B 
S 


i _ 11 - + - cos cos (co - g)t + 1 Lz°2*4L 


_ .T* -r?0. 

Va 

-2 (b S °) 

-2 (B ■ S ) 

V x xj 

S 2 


S 2 


2 • fi)t + — - — ~cos (co + fi)t 


V 2 


Expanding these terms to second order in L/fi and (co - fi)t, we find again 
that B/S « 1 and 


B 

S 


0 = “ = [L 2 t] 2 + 


(i 1 )' 


f - cos (co + fi)t 
2 


■Va 


( 20 ) 


and again we have an oscillatory term in the expression for © but the detailed 


time dependence is different from case i. For t » (77) > ® again becomes 

linear in t. 2 


© 


3. S^ 0 = S, S x ° = 0, S ° = 0; i. e. , S is initially parallel to the 

earth's axis. Equations (10) give ip = 0, Cj = S cos <p, C 2 = S sin <p, 


S = -S 
x 


sin (p cos cp sin fit + sin <p + ° os ft) s ^ n 


. , (1 - cos <p) . , „ VJ . 

+ sm <p — — — ^ sm (co + fi)t 

z 


10 



y 


• . , ^ , (1 + cos (b) 

sin <p cos <p cos fit - sm <p r - / - cos (co - fi)t 

u 


, (1 - cos <b) 

+ sm (p — — r — cos (co + fi) t 


S = S [cos 2 <b + sin 2 <±> cos cot] 
z 


The dominant terms in S become, to first order in L/fi, 


S = -SI— 1 ) sin fit 
x \fi 


( 21 ) 


S = Sl-^ ) [1 - cos fit] 
y Vfi 


( 22 ) 


with only terms with frequency fi. S is constant to first order in L/fi. The 

z 

expression for the precession angle is 


e « B - 
0 s 




y 2 


i/ 

= \T2 sin 0 [ 1 - cos cot] 1 


( 23 ) 


so that 0 is periodic with frequency co ~ fi . 


THE MOTIONAL PRECESSION 


The complete version of Schiff's equation, including both geodetic and 
motional precession, is 


dS — — 

— = S x L' 
dt 


(24) 


v, _ r gi 

L' — L + 9 — r 

c 2 r J 


3r \ — 

(v ' r) - (7 


L + H 
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where a is the earth's angular velocity and I its moment of inertia. Now the 
first term in H oscillates with period equal to the orbital period so that the 
exact solution to equation (24) for a precessing circular orbit is much more 
complicated than the solution to equations (5), which we have already found. 
However, the effect of H can be taken into account by replacing H by its 
average over one orbital period T. The justification for this statement is 
contained in the following argument: 1 If we integrate equation (24) over an 
orbital period t and divide by T, we obtain 


S(t + T) - S(t) 
T 


t+T 

-J dt'[S(t') x L'(t')] 


Now S (t 1 ) can be expanded in a Taylor series about t: 

— — ds (t) 

S(t') = S(t) + — (V - t) + . . . 

= S(t) + [S(t) x L'(t)J(f - t) + ( ter “ s of . M £ er 

I order m L' 

Defining the average value of S over an orbital period as 

x v i t+T 

= - f dt* S(t') , 

t 


(25) 


we see that if we differentiate this equation we obtain the left hand side of 
equation (25). So that substituting the Taylor series for S (t 1 ) in the right 
hand side of equation (25) gives 

d<^S(t)> 1 t+T 

dt = “t / dt ' L' (t') x {lf(t) + [S(t) x L'(t)] (t' - t) + . . . 

t 


Now if we neglect the second and higher terms in the bracket and we replace 
S(t) by <S(tj)> in the bracket, this equation will still be valid to first order 
in L' , and it can then, be written 


d 




1 . A more rigorous argument could probably be given by applying some kind 
of classical perturbation theory to equation (24), 



where 


t+T _ 

<L'(t)> = - / dt'L' (f) . 

t 

This means that to first order in L' (t) , the average value of S (t) over an 
orbital period satisfies equation (24) with L' (t) replaced by <^L I (t)^> so 
that we can replace H by in equation (24) if we are interested only in 

the average value of S over an orbital period. 

This type of argument is borne out by the solution given in the first 
section if we consider T to be t, the period of the precession of L rather than 
the orbital period. The terms with period r = 2-k/Q, in equations (5) are seen 
to contribute nothing to the average value of S, which is seen to precess about 
the earth's axis; i. e. , about the average value of L, to first order in L. 


We also note that this argument allows us to drop the restriction to 
circular orbits in the previous calculation. Since r” x "v is a constant for 

noncircular orbits, we simply replace pj- by , the average of p$- over 

the orbital period, in the expression for L. This is given by [2, equation 
11.15.4] 



where a is the semimajor axis of the orbit and e the eccentricity. 


We will now compute <^H)> for a circular orbit with inclination 

angle i. To do this, we construct a coordinate system such that the orbit 
is the x'-y' plane and the orbit normal is the z' axis. We will let"a be in 
the y'-z' plane, so that it will make an angle i with the z' axis. In this 
coordinate system the orbit is given by 

r ' = -a sin at r ' = 0 

x z 


r ' = a cos at T = — 

y a 


and a ' = 0, a ' = cr sin i, a ' = a cos i, so that (a • r) = acr ’ cos at. 
x y z y 
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Letting <^H J^> , </h^> , </H ; ^> deaote the components of < \H/ > 

in this coordinate system we have 


< H x'> 


GI i 


T 


cV T 


I [ (-3 sin at) (a 1 cos at) ]dt 

“ y 


< h ;> 


-3GIa 1 

y 



27T 

sin 2 a t 

a 

2a 



o 


= 0 


GI 1 

f 1(3 cos at) (a ' cos at) 


a »]dt 

y 


<-„•> - 


GI , 


~ 2 "i a 1 

c^ar y 

T 

GI 

_ f 

2c 2 a 3 y 

GI 1 

1 

r 

(?£? T 

J 


3 l at 1 . _ A 

— ( — + 7 sm 2a t ) 
a y 3 4 ) 


2n 

a 


- T 


-IG , 
H cr 1 
c^a 15 z 


Performing a rotation about the x'-axis through an angle i we have 
< H y) = <»,’> ° 0S ‘ - < H z> Bln 1 


3GIcr r . . , 

2^r I sin 1 cos ^ 


< H z> = < H y’> Sini+ < H z’> cos i 


Glcr 

c 2 a^ 


sm" i 2 • 

— - — - COS 1 


5S? [1 - 3cos2il 


< H X> 
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Since the orbital angular momentum, and thus L, remains in the y-z plane after 

the rotation, we see that has components only in the plane which contains 

L and a. The components of <^H )> perpendicular to this plane average to zero, 

and we expect this to remain true for noncircular orbits. Now as the orbit 

plane and L precess, *\H will remain in the plane described by L and’a, 

and L + <^H /■ will precess about the earth axis with period t = This 

means that the solution to equation (24) for L' = L + { H ^ can be obtained by 

substituting L' for L in the solution which we found previously for equation (1) . 
In other words we simply substitute Lj' for and L 2 ' for L 2 in equations (5) 
and all subsequent equations in the first section where 

T , t 3GLt r -i 

Lj 1 = Lj + 2 .3 l sln 1 cos G 

£ c r 

W = L 2 + I 1 - 3 cos^ x] . 

It is now cleaxythat H contributes to the precession no matter what 
initial orientation for S is chosen. If S is initially oriented pex-pendicular to 
the earth's axis (cases 1 and 2) , it will, apart from effects with period t, 
precess about the z-axis with angular velocity L 2 * as described by equations 
(16), (17), (18), (19), and (20) with Lj (L 2 ) replaced by Lj 1 (L 2 '). 

Both L and H will contribute to the precession. If S is initially 
aligned parallel to the earth's axis (case 3) , it will exhibit oscillatory 
behavior with period r, _as described by equations (21) , (22) , and (23) , with 
Lj replaced by L^. So H contributes to the precession in this case also. 

For a general inclination angle i, H will contribute to the precession for any 
initial orientation of S, since both L^' and L 2 ' contain the effects of H. (The 
only exceptions to this statement occur for a polar ox-bit, an equatorial orbit, 
and an orbit for which cos i = 1A/TT, i = 54 degrees. ) So that in general it 
is not possible to separate the effects of L and H and thus verify the existence 
of the geodetic and Lense-Thirring precessions separately for an inclined orbit. 

SKYLAB II ORBIT 


The nominal Skylab II orbit will be^circular at 235 nautical miles 
inclined at 35 degrees. We will compute S for this case for the three initial 

orientations discussed pi-eviously. For this orbit 

r = 6. 8 x 10® meters 
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and Kepler's Law gives for the period r 


T = 2it 



MG 


93 min 


= (27r)/T 


We have R/r = 0. 94, so L = 7.2 arc sec/year, and L 2 = L cos i = 5. 9 arc 
sec/year. 

Now in units of earth radii R 
J 2 = 1.08 x 10" 3 

so that equation (3) gives 0 = 6. 5 degrees/ day, t = 55 days, and 
^ sin i = 0. 10 arc sec. 

The components of S for cases 1, 2, and 3 are plotted in Figures 1 
through 7, using equations (16), (17), (18), (19), (20), (21), (22), and 
(23). For this case, L 2 differs from L 2 ’ by less than 1 percent. 




Figure 1. Case 1: — , 0 
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